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Realce 

geométrico 


Magnificación 
con  interpolación 


Operaciones  de  magnificación  de  la  imagen. 


Introducción 


Se  entiende  como  realce  geométrico  un  conjunto  de 
operaciones  destinadas  a  mejorar  la  calidad  visual  de 
la  imagen  alterando  algunas  de  sus  características 
geométricas,  sin  alterar  sustancialmente  a  su  radio¬ 
metría.  Este  objetivo  puede  abordarse  por  diversos 
caminos: 

-i  Modificando  el  tamaño  de  los  objetos  en  la  salida 
gráfica  de  la  imagen. 

-i  Minimizando  defectos  localizados  o  generalizados 
de  la  imagen. 

-i  Realzando  rasgos  que  mejoren  la  percepción  de 
los  elementos  de  la  imagen. 


Estas  operaciones  de  realce  geométrico  son  útiles 
para  el  análisis  visual  de  la  imagen  sobre  pantalla  o 
sobre  papel.  Son  deseables  e  incluso  necesarias  en 
numerosas  ocasiones  de  cara  a  la  producción  carto¬ 
gráfica  o  a  la  presentación  de  los  resultados  de  ope¬ 
raciones  realizadas  con  imágenes  sobre  cualquier 
tipo  de  soporte,  especialmente  las  dos  últimas.  Sin 
embargo,  no  deben  llevarse  a  cabo  si  tras  ellas  se  pre¬ 
tende  obtener  resultados  con  valor  físico,  puesto  que 
comportan  la  alteración,  siquiera  mínimamente,  de 
los  ND  originales. 

Entre  las  operaciones  de  realce  geométrico,  una  de 
las  más  empleadas,  aunque  no  la  única,  es  sin  duda 
el  filtrado  digital. 


Cambios  de  escala 


Es  impropio  aplicar  el  concepto  de  escala  a  la  imagen 
digital,  porque  sólo  es  un  registro  electrónico  alma¬ 
cenado  en  forma  de  archivo,  incluso  sin  estructura 
matricial  sino  lineal.  La  imagen  digital,  pues  no  tiene 
escala.  Sólo  tiene  sentido  hablar  de  escala  cuando  esa 
imagen  es  representada  en  algún  soporte  plano  e  in¬ 
deformable. 

En  general  se  prefiere  utilizar  el  concepto  de  resolu¬ 
ción,  sobreentendida  ésta  como  espacial,  es  decir,  el 
tamaño  nominal  de  la  celda  en  el  terreno  y  en  el  nadir 
de  la  imagen.  Corresponde  a  las  denominaciones  de 
GSD  ( Ground  Sample  Distance,  distancia  de  mues- 


treo  en  tierra),  GIFOV  ( Ground-projected  Instanta- 
neous  Field  of  View,  campo  de  visión  instantáneo 
proyectado  en  tierra).  Este  texto  utilizará  preferente¬ 
mente  la  denominación  de  tamaño  de  la  celda,  enten¬ 
diendo  ésta  en  el  terreno. 

Existen  dos  tipos  de  cambio  de  escala: 

-i  La  magnificación  o  zoom.  Consiste  en  replicar 
cada  celda  una  o  más  veces  de  modo  que,  del  pixel 
original  de  pantalla  que  representa  cada  celda  en 
una  presentación  natural  de  la  imagen,  cada  celda 
pase  a  estar  representada  por  n2  píxeles,  siendo  n 
el  grado  de  magnificación.  Es  importante  insistir 
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en  que,  aunque  la  celda  se  perciba  con  mayor  ta¬ 
maño  en  la  pantalla,  la  resolución  espacial  sigue 
siendo  la  misma.  No  se  incrementa  el  nivel  de  de¬ 
talle,  solamente  se  han  repetido  valores  de  celda  y 
se  han  utilizado  píxeles  contiguos  para  simular 
una  celda  más  grande.  Ciertos  sistemas  de  trata¬ 
miento  digital  de  imágenes  disponen  de  una  mo¬ 
dalidad  de  magnificación  con  interpolación  que 
permite  suavizar  las  transiciones  entre  los  valores 
replicados  de  modo  que  se  reducen  los  saltos  brus¬ 
cos  de  valores.  Salvo  que  la  variable  representada 


El  filtrado  digital 

El  filtrado  digital  es  una  transformación  de  la  imagen 
mediante  operadores  locales  que  involucran  a  la 
celda  que  se  pretende  modificar  y  a  un  entorno  cer¬ 
cano  de  ella.  El  nuevo  ND  de  la  celda  filtrada  será 
función  no  sólo  del  valor  de  la  celda  original,  sino 
también  de  los  ND  de  las  celdas  contiguas. 

Los  filtros  son  operadores  que  se  aplican  con  objeti¬ 
vos  opuestos,  dando  origen  a  dos  grandes  familias, 
los  filtros  de  paso  bajo,  que  suavizan  la  imagen  final, 
aproximando  los  ND  originales  a  los  de  sus  vecinos, 
bien  para  destacar  elementos  de  interés,  y  los  filtros 
de  paso  alto,  que  exageran  las  diferencias  y  aumentan 
el  contraste  local. 

Los  filtros  digitales  pueden  ser  lineales  y  no  lineales. 

El  filtrado  digital  es  una  operación  de  convolución. 
Se  define  la  convolución  de  dos  funciones  unidimen¬ 
sionales  f\(t)  yfzit )  como  el  resultado  de  integrar  el 
producto  de  las  dos  funciones  después  de  desplazar 
una  de  ellas  una  distancia  x,  que  varía  entre  -oc  e  ce, 
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en  la  imagen  tenga  alta  autocorrelación,  el  resul¬ 
tado  obtenido  no  tiene  por  qué  ser  más  fiel  a  la 
imagen  original. 

-i  La  reducción.  Consiste  en  prescindir  de  una  o  más 
celdas  siguiendo  un  patrón  de  elección,  de  modo 
que  resulte  una  imagen  reducida  compuesta  por 
una  muestra  de  la  original.  Esta  operación  es  ne¬ 
cesaria  para  poder  tener  una  visión  de  conjunto 
de  una  imagen  cuya  extensión  real  en  número  de 
celdas  (columnas  x  líneas)  supere  el  tamaño  del 
monitor  o  de  la  salida  gráfica  que  se  desea  obte¬ 
ner. 


es  decir: 


=  )  ®fi(t)=  J  fiOf2{t~x)dx  (1) 


en  donde  0  es  el  símbolo  de  la  convolución. 

El  filtro  digital  es  una  operación  de  convolución  en¬ 
tre  la  imagen  original  z(x,y)  y  la  función  filtro  h(x,y). 
Pero  para  que  esta  operación  sea  aplicable  dentro  del 
ámbito  de  las  imágenes  digitales  se  requieren  tres  re¬ 
quisitos: 

-i  Debe  aplicarse  bidimensionalmente. 

—i  Debe  ser  una  función  finita. 

-i  Debe  ser  una  función  discreta. 

La  convolución  (1)  se  puede  enunciar  de  modo  bidi- 
mensional,  que  es  más  apropiada  para  el  ámbito  del 
tratamiento  de  imágenes.  La  convolución  entre  dos 
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funciones  bidimensionales  h(x,y )  y  z(x,y)  es: 
z'(x,y)  =  h(x,y)®z(x,y)  = 

+00 

=  ||  h(x,  y)f[x-u,  y-v) dudv 


El  filtrado  digital. 


donde  la  variable  auxiliar  móvil  x  de  la  ecuación  (1) 
se  ha  cambiado  por  las  variables  u  y  v. 

Para  cumplir  la  segunda  y  tercera  condición  es  nece¬ 
sario  cambiar  los  límites  de  integración  y  acotarlos 
según  la  dimensión  de  la  imagen  y,  además,  ya  que  la 
imagen  digital  es  discreta  y  no  continua,  es  preciso 
sustituir  la  doble  integral  por  una  doble  sumatoria. 
Normalmente,  el  filtro  consiste  en  una  matriz  desli¬ 
zante  de  dimensión  impar  ( 2g+\ )  x  (2/z+l)  y  su  con- 
volución  con  la  imagen  se  expresa: 

z'(m,n )  =  w(&,/)®  f  (m,n)  = 

s  *  (3) 

=llHw(k’l)f(m+k  ,n+l ) 


Donde  donde  m  y  n  son,  respectivamente,  las  colum¬ 
nas  y  líneas  de  la  imagen,  y  k  y  /  son  las  columnas  y 
líneas  de  la  matriz  de  filtrado. 


wu 

wn 

wn 

w21 

W22 

W22 

W31 

W22 

w33 

sobre  una  imagen  z(m,n)  da  lugar  a  la  imagen: 

z'(m,n)  =  wnz(m-\,n-\)  +  wX2z{rn,n-\)  +  wuz{m+\,n-\)  + 

+  w21z(m- \,n)  +w22z{m,n)  +  w23z(m+l,«)  + 

+  w3lz  (m-1,  n+ 1)  +  w32z  (m,  n+ 1)  +  w33z  (m+1,  n+ 1) 

Una  vez  obtenido  el  valor  de  z\m,rí),  se  repite  la  ope¬ 
ración  para  las  posiciones  ( m+\,n ),  ( m+2,n ),  etc. 
hasta  terminar  con  la  línea  n  y  seguir  con  la  n+1,  y 
así  sucesivamente  hasta  recorrer  la  imagen  completa. 

Conviene  que  la  matriz  de  filtrado  sea  normalizada 
para  no  variar  la  luminancia  general  de  la  imagen  de 
salida.  La  normalización  consiste  en  hacer  que  la 
suma  de  todos  sus  elementos  sea  igual  a  la  unidad. 
Si  la  suma  de  los  elementos  de  peso  de  la  matriz  de 
filtrado  no  cumple  esa  condición,  se  normaliza  divi¬ 
diendo  el  resultado  por: 

HWU  (5) 


Los  elementos  w  de  la  matriz  de  filtrado  se  denomi¬ 
nan  coeficientes  de  peso.  El  entorno  [-g,  g]  •  [- h ,  h]  es 
la  ventana  del  filtro,  en  la  cual  usualmente  g  =  h,  es 
decir,  la  ventana  suele  ser  cuadrada.  Según  sea  el  va¬ 
lor  de  g  y  h,  la  influencia  del  entorno  de  las  celdas  en 
el  ND  final  podrá  ser  de  mayor  o  menor  alcance.  Los 
filtros  más  utilizados  son  de  3x3  elementos,  y  en  al¬ 
gunos  casos  de  5x5,  aunque  pueden  aplicarse  filtros 
mayores  en  determinados  casos. 

La  aplicación  de  un  filtro: 


El  filtro  puede  ser  definido  por  el  usuario,  de  acuerdo 
a  sus  necesidades.  Según  la  naturaleza  y  el  efecto  que 
produzcan  sobre  la  imagen,  los  filtros  se  clasifican  de 
diversas  maneras  Entre  los  filtros  lineales  se  encuen¬ 
tran: 


-i  Filtros  de  paso  bajo, 
-i  Filtros  de  paso  alto. 
— i  Filtros  de  muesca. 

-i  Filtros  direccionales. 
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— i  etc. 

Entre  los  no  lineales,  los  más  habituales  son: 

-i  Filtros  de  moda  o  de  mayoría. 

-i  Filtros  de  mediana. 

—i  Filtros  de  refuerzo  o  de  extracción  de  bordes, 
-i  etc. 

Filtros  de  paso  bajo 

Los  filtros  de  paso  bajo  son  filtros  que  suavizan  la 
imagen  uniformizando  cada  ND  original  en  relación 
con  los  de  su  entorno.  La  denominación  de  paso  bajo 
es  una  reminiscencia  del  tratamiento  de  señales  elec¬ 
trónicas  unidimensionales,  del  cual  deriva  el  trata¬ 
miento  digital  de  imágenes  habitual  en  teledetección. 
Las  señales  electrónicas  son  manipuladas  en  el  domi¬ 
nio  de  la  frecuencia  y,  dentro  de  él,  los  filtros  se  des¬ 
tinan  a  eliminar  los  componentes  de  altas  o  de  bajas 
frecuencias,  que  en  correspondencia  con  los  sonidos 
transformados  en  señal  equivalen  a  sonidos  agudos  y 
graves  respectivamente.  Los  filtros  de  paso  bajo  eli¬ 
minan  las  altas  frecuencias  y  dejan  pasar  solamente 
las  bajas  y  medias. 

Ejemplos  de  filtros  de  paso  bajo  son  los  siguientes: 


'1  1  f 

1  1  1 

1 

'i  i  f 

1  2  1 

1 

1 

LJ  H- 

4^  tO 

LJ  ^ 
_ l 

1  1  1 

7o 

1  1  1 

16 

1  2  1 

(6) 

El  primero  de  ellos  se  denomina  filtro  de  media  por¬ 
que  devuelve  el  valor  medio  calculado  entre  el  cen¬ 
tral  y  los  ocho  más  próximos.  Los  otros  dos  filtros 
atribuyen  más  importancia  al  ND  central  que  a  su  en¬ 
torno  inmediato  y,  por  tanto,  en  la  imagen  de  salida 


se  percibirá  un  efecto  de  suavizado  menos  intenso 
que  con  el  filtro  de  media. 

Se  aplican  para  eliminar  ruido  de  la  imagen  y  para 
suprimir  el  bandeado.  Se  entiende  por  ruido  el  con¬ 
junto  de  errores  aleatorios  introducidos  en  la  imagen 
durante  el  proceso  de  adquisición,  transmisión  o  gra¬ 
bación  y  puede  materializarse  en  píxeles  anormal¬ 
mente  claros  u  oscuros  salpicados  de  una  forma  más 
o  menos  aleatoria,  efecto  que  se  conoce  con  el  nom¬ 
bre  de  sal  y  pimienta. 

Este  tipo  de  filtros  son  también  utilizados  para  suavi¬ 
zar  contrastes  entre  diferentes  zonas  de  la  imagen  y 
para  reducir  la  variabilidad  de  la  imagen  antes  de  pro¬ 
ceder  a  una  clasificación  digital.  Presentan  en  gene¬ 
ral  el  inconveniente  de  difuminar  excesivamente  los 
bordes. 

Filtros  de  paso  alto. 

Son  operadores  que  realzan  los  rasgos  lineales  de  la 
imagen  aumentando  los  contrastes  locales.  Se  llaman 
así  porque  refuerzan  los  componentes  de  alta  fre¬ 
cuencia  en  la  imagen,  aunque  pueden  no  eliminar  los 
de  baja  frecuencia. 

El  procedimiento  es  el  contrarío  al  seguido  anterior¬ 
mente:  se  trata  ahora  de  aumentar  digitalmente  las 
diferencias  entre  ND  vecinos.  Existen  dos  formas  de 
obtener  un  refuerzo  de  las  zonas  de  alta  variabilidad: 
un  procedimiento  sencillo,  aunque  poco  utilizado 
consiste  en  sumar  a  la  imagen  primitiva  la  diferencia 
entre  ella  misma  y  la  producida  por  un  filtro  de  paso 
bajo: 

kj  ~  kj  +  ( kj  ~  kj  )  —  —  Z¡j  ^ 

donde  z”  es  la  imagen  realzada,  z  la  original  y  z’  es 


La  imagen  central  es  la  parte  del  centro  del  rose¬ 
tón  superior.  La  imagen  inferior  es  el  resultado 
de  aplicarle  un  filtro  de  media. 
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Resultado  de  aplicar  un  filtro  de  paso  alto. 


la  imagen  filtrada  con  un  filtro  de  paso  bajo. 

La  otra  forma  es  aplicar  filtros  de  convolución  de  re¬ 
fuerzo  de  bordes.  Para  realzar  las  fronteras  entre 
áreas  contiguas  de  la  imagen  independientemente  de 
su  dirección,  se  necesita  que  el  filtro  tenga  un  gra¬ 
diente  isótropo  de  pesos  en  las  ocho  direcciones  po¬ 
sibles  a  partir  del  elemento  central  de  la  matriz,  esto 
es,  que  los  ocho  elementos  adyacentes  al  central  ten¬ 
gan  pesos  idénticos.  Estos  filtros  se  caracterizan  por 
contener  valores  negativos  y  positivos.  Si  solamente 
se  desea  detectar  los  bordes  sin  mantener  la  imagen, 
la  suma  de  todos  los  coeficientes  periféricos,  debe  ser 
igual,  pero  de  signo  contrario  al  peso  asignado  a  la 
celda  central. 

Algunos  filtros  comúnmente  usados  son: 


“  0 

-1 

0" 

“  0 

-1 

0" 

-1 

4 

-1 

-1 

5 

-1 

0 

-1 

0 

0 

-1 

0 

de  los  cuales,  el  primero  es  conocido  como  filtro 
laplaciano.  Obsérvese  que  la  suma  de  sus  elementos 
de  peso  es  cero.  Esto  quiere  decir  que  el  valor  del 
brillo  final  de  la  imagen  filtrada  será  muy  reducido. 
La  imagen  habrá  desaparecido  y  en  su  lugar  sola¬ 
mente  destacarán  los  bordes  que  en  ella  hubiera.  Si 
se  desea  conservar  la  imagen,  pero  enriquecerla  con 
los  rasgos  lineales  realzados  se  puede  sumar  la  ima¬ 
gen  original  al  resultado  del  filtro  laplaciano.  Pero 
como  el  filtrado  es  una  operación  lineal,  lo  mismo  da 
sumar  las  imágenes  al  final  que  sumar  al  filtro  la  fun¬ 
ción  identidad,  que  sería: 


0  0  0 
0  1  0 
0  0  0 


(9) 


Su  suma  con  el  laplaciano  origina  el  segundo  de  los 
filtros  (8). 

Otros  posibles  filtros  pueden  ser: 


“-1  -1  -f 

“1-2  f 

“-1  0  -f 

-1  9  -1 

-2  5  -2 

0  5  0 

-1  -1  -1 

1  -2  1 

-1  0  -1 

(10) 

cuya  selección  dependerá  del  valor  que  pretenda  dár¬ 
sele  al  ND  original  de  cada  posición  en  relación  con 
los  circundantes. 


Filtros  direccionales 

Los  filtros  direccionales  son  operadores  de  gradiente 
anisótropo  que  tienden  a  realzar  los  rasgos  lineales 
de  la  imagen  primando  determinadas  direcciones.  El 
efecto  de  refuerzo  lineal  puede  ser  unidimensional: 

z'(m,n)  =  z(m,n)~  z^1  (m,n)  (11) 

cuya  expresión  significa  sustraer  de  sí  misma  la  ima¬ 
gen  desplazada  una  posición  en  cualquier  dirección 
elegida.  Si  no  hay  cambios  en  esa  dirección,  el  resul¬ 
tado  para  la  celda  será  cero,  y  será  tanto  más  alto  en 
valor  absoluto  cuando  mayor  diferencia  de  ND  haya 
en  dicha  dirección. 

Este  reforzador  lineal  también  puede  operar  sobre 
los  dos  sentidos  de  la  dirección  elegida: 


Imagen  real  y  filtro  laplaciano. 
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z'(m,n)  =  [z(m,n)- z'*1  (m,n)  +  z(m,n)~  z^  1  (m,n) 


(12) 


Igual  que  antes,  valores  de  z '  próximos  a  cero  indi¬ 
can  zonas  de  escasa  variación  y  viceversa,  valores  al¬ 
tos  indican  que  ha  habido  un  cambio  local  de  lumi- 
nancia,  es  decir,  que  existe  un  borde. 


Filtros  no  lineales 

Existen  filtros  que  no  mantienen  la  propiedad  de  li- 
nealidad  matemática,  fundamentalmente  porque  sus 
elementos  de  peso  son  variables  con  la  posición  sobre 
la  que  se  aplican.  En  estas  circunstancias,  aventurar 
el  resultado  del  filtrado  puede  ser  difícil  sin  probar 
antes. 


Sin  embargo,  el  filtrado  direccional  más  utilizado  está 
basado  en  la  convolución  mediante  los  filtros  de  Pre- 
witt: 


'1 

1 

f 

1 

-2 

-1 

1 

-1 

-1 

NO 

"1 

1 

-f 

1 

-2 

-1 

1 

1 

-1 

O 

"1 

-1 

-f 

1 

-2 

-1 

1 

1 

1 

SO  s 


1 

1 

f 

-1 

-2 

1 

-1 

-1 

1 

NE 

'-1 

1 

f 

-1 

-2 

1 

-1 

1 

1 

E 

'-1 

-1 

f 

-1 

-2 

1 

1 

1 

1 

SE 


(13) 


Uno  muy  difundido  es  el  filtro  de  mediana.  Se  co¬ 
noce  también  como  ventana  de  Tukey,  y  fue  ideado 
para  suprimir  ruido  de  la  escena  sin  difuminarla  de¬ 
masiado  y  respetando  así  los  bordes.  Consiste  en  asig¬ 
nar  a  la  posición  central  de  la  ventana  de  filtrado,  que 
suele  ser  de  3x3,  un  valor  igual  a  la  mediana  estadís¬ 
tica  de  los  ND  incluidos  en  ella.  No  existe,  pues,  una 
ventana  de  filtrado  única  para  toda  la  imagen,  sino 
que  ésta  cambia  de  coeficientes  según  se  desplaza  a 
través  de  la  escena.  Para  una  ventana  de  3x3,  el  filtro 
sería: 


z(m-\,n—\)  z(m-l,n-l)  z(m-l,n-l) 
z(m-l,n)  z(m,n )  z(m+\,n) 

z(m-l,n+l)  z(m,n+ 1)  z(m+l,«+l) 

(14) 

Tiene  un  efecto  de  paso  bajo  y  normalmente  respeta 
mejor  las  formas  que  el  filtro  de  media,  pero  al  tra¬ 
tarse  de  un  operador  no  lineal,  sus  resultados  no  son 
previsibles  de  forma  general. 


z!{m,n )  =  Med 


Que  realzan  los  bordes  detectados  en  las  direcciones 
norte,  sur  este,  oeste  e  intermedias,  dependiendo  del 
gradiente  de  sus  elementos.  Obsérvese  que  siempre 
guardan  una  simetría  respecto  a  un  eje  perpendicular 
a  la  dirección  en  la  que  detectan  el  borde. 


Otro  muy  utilizado  es  el  filtro  de  moda  o  de  mayoría. 
A  la  posición  central  de  la  ventana  de  filtrado  se  le 
asigna  un  valor  igual  a  la  moda  estadística  de  los  ND 
incluidos  en  ella.  No  existe  tampoco  en  este  caso  una 
ventana  de  filtrado  única  para  toda  la  imagen,  sino 
que  ésta  cambia  de  coeficientes  según  se  desplaza  a 
través  de  la  escena.  Para  una  ventana  de  3x3,  este 


Resultados  de  aplicar  filtros  de  Prewitt  N,  S,  E  y  O. 
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Imagen  original  y  filtro  de  mediana. 


filtro  puede  enunciarse  del  siguiente  modo: 

z{jn-\,n—\)  z(m-l,n-l)  z(m-l,n-l) 
z(m-l,n)  z(m,n )  z(m+l,«) 

z(m-l,n+l)  z(m,n+ 1)  z(wi+l,rc+l) 

(15) 

Su  efecto  es  el  de  eliminar  píxeles  sueltos  de  un  valor 


z'(m,n )  =  Mod 


diferente  al  de  su  entorno.  Se  aplica  fundamental¬ 
mente  sobre  imágenes  cualitativas.  Es  frecuente  uti¬ 
lizarlo  sobre  el  resultado  de  una  clasificación  digital 
para  obtener  un  aspecto  más  generalizado  de  la  dis¬ 
tribución  de  clases,  ofreciendo  grandes  zonas  homo¬ 
géneas  que  el  usuario  puede  interpretar  más  fácil¬ 
mente,  aun  a  riesgo  de  perder  una  cuota  de  veracidad 
en  los  resultados. 

Existe  un  gran  grupo  de  filtros  no  lineales  destinados 
a  detectar  contornos.  Para  reconocer  automática¬ 
mente  los  contornos,  estos  operadores  calculan  las 
diferencias  entre  los  ND  de  unja  celda  y  las  contiguas 
en  todas  las  direcciones.  El  trazado  de  los  contornos 
se  basa  en  representar  los  vectores  direccionales  cu¬ 
yos  módulos  son  las  diferencias  en  todas  las  direccio¬ 
nes  que  superen  un  determinado  umbral  seleccio¬ 
nado.  Precisamente  en  la  elección  de  los  umbrales 
depende  en  gran  medida  la  eficacia  de  estos  filtros. 
Pueden  aparecer  bordes  no  deseados  o  perderse 
otros  realmente  buscados. 

El  filtro  de  Roberts  se  trata  de  un  filtro  no  lineal  apli¬ 
cado  mediante  una  ventana  de  2x2: 

z'  {i,  j )  =  ^[z  (i,  j )  -  z  (i+1,  ;'+l)]2  +  [z  (i,  ;'+l)  -  z  (i+ 1,  ;)]2 

(16) 

o  bien: 

z(i,  j )  =  |  z(i,  j)~z(i+ 1,  j+l)  +  z(i,  7+1)  -  7 (i+1,  j)  | 

(17) 


Resultados  de  aplicar  un  filtro  de  moda  sobre  la  imagen 
superior 
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Filtrado  digital 


Para  la  extracción  de  bordes  también  se  aplica  el  fil¬ 
tro  de  Sobel: 

z'{i,j)  =  4cr+72  (18) 

siendo  la  matriz  representada  una  ventana  de  la  ima¬ 
gen  de  3x3  centrada  sobre  la  celda  (i,  j )  y  donde  los 
coeficientes  C  y  F  se  calculan  del  siguiente  modo: 


C  =  (z2+2z3  +  z4)-(zo  +  2z7  +  z6) 

F  =  (zo  +  2zi+z2)-(z6+2z5  +  z4) 

(19) 

Para  usar  el  filtro  de  Kirsh  se  definen  las  variables 
variables  Sn  y  Tn: 

Sn=Zn+  4+1  +  4+2  ^Q) 

T„  =  4+3  +  4+4  +  4+5  +  4+6  +  4+7 

siendo  n  cualquiera  de  las  ocho  posiciones  señaladas 
alrededor  de  la  celda  (ij): 


z0  zl  z2 
4  z[i,j]  z3 
z6  z5  z4 


Zl  ^2  ^3 

4  z[i,j]  4  (21) 

z7  z6  z5 

Se  calculan  los  ocho  valores  de  Sn  y  Tn  para  n  va¬ 

riando  entre  0  y  7  y  se  adopta  el  valor  del  filtro  como: 

AT(f,2)  =  M«r[|5S.-3r„|]  (22) 

Finalmente,  otro  de  los  filtros  de  extracción  de  bor¬ 
des  utilizado  con  relativa  frecuencia  es  el  de  Wallis. 
Este  operador  considera  una  celda  de  imagen  como 
borde  si  el  logaritmo  de  su  ND  excede  de  la  media  de 
los  logaritmos  de  los  cuatro  vecinos.  Asumiendo  la 
nomenclatura  de  (21): 


Original  y  filtro  de  Roberts. 


Resultado  del  filtro  de  Sobel. 


Original  y  filtro  de  Kirsh. 
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z'(U)  = 

log[z(¿,;)]-j[log(z2)  +  log(z4)  +  log(z6)  +  log(z8)]  = 

log  [  z  ( /,  j )]  -  i  l°g  ( z2  z4  76  78 ) 

o,  lo  que  es  lo  mismo: 


ÍMutfl 


(23) 


Transformada  de  Fourier  de  una  imagen  mos¬ 
trando  las  frecuencias  espaciales  en  una  escala 
de  0  a  1. 


Filtrado  en  el  dominio  de  la  frecuencia 


Cualquier  señal  periódica  puede  enunciarse  como 
una  sumatoria  infinita  de  funciones  periódicas  de  dis¬ 
tintas  frecuencias.  Este  nuevo  espacio  de  referencia 
para  la  descripción  de  la  función  se  conoce  como  do¬ 
minio  de  la  frecuencia.  La  serie  de  Fourier  permite 
realizar  el  cambio  de  dominio  temporal  al  de  la  fre¬ 
cuencia.  Una  función  periódica  en  forma  de  serie  de 
Fourier  queda  unívocamente  definida  mediante  la  es¬ 
pecificación  de  los  coeficientes  de  dicha  serie.  La  re¬ 
presentación  gráfica  de  los  coeficientes  complejos  de 
la  serie  en  función  de  la  frecuencia  angular  se  deno¬ 
mina  espectro  de  amplitud  y  la  representación  gráfica 
del  ángulo  de  fase  en  función  de  la  frecuencia  es  el 
espectro  de  fase: 

La  generalización  del  análisis  de  Fourier  y  su  utiliza¬ 
ción  en  el  campo  de  la  teledetección  requiere  desa¬ 
rrollar  un  método  aplicable  a  funciones  no  periódi¬ 
cas,  a  las  cuales  puedan  ser  asimiladas  las  imágenes 
digitales.  La  integral  o  transformada  de  Fourier  per¬ 
mite  describir  funciones  no  periódicas  en  el  dominio 
de  la  frecuencia.  De  forma  análoga,  la  transformada 
inversa  de  Fourier  permite  pasar  de  nuevo  al  dominio 
temporal  original. 

En  el  análisis  digital  de  imágenes  todavía  hay  que 


realizar  algunas  adaptaciones  más  para  poder  aplicar 
la  transformada  de  Fourier.  En  primer  lugar,  las  fun¬ 
ciones  deben  ser  finitas  y  discretas,  como  lo  son  las 
imágenes  digitales.  Además,  la  variable  temporal 
debe  ser  reemplazada  por  la  espacial,  como  corres¬ 
ponde  a  las  imágenes  de  satélite,  en  contraste  con  las 
señales  electrónicas  o  de  audio.  Por  último,  tampoco 
son  útiles  las  funciones  unidimensionales,  puesto  que 
la  imagen  se  desarrolla  en  el  plano  bidimensional. 
Por  lo  tanto,  debe  reconsiderarse  la  teoría  y  adaptar 
el  análisis  de  Fourier  al  caso  de  las  imágenes  digita¬ 
les.  Para  ello  se  enuncia  la  transformada  discreta  y 
bidimensional  de  Fourier. 

Sea  (|)(i,j);  i,  j  =  0,  -,  K- 1,  el  brillo  de  la  celda  (i,j)  de 
una  imagen  cuadrada  de  K2  filas  por  columnas. 

Su  transformada  discreta  de  Fourier  es: 

«M-ZZ+OJ)  <24> 

i=0  j=0 

En  esta  expresión  i  =  4-1  y  ry  s  son  las  frecuencias 
espaciales  en  las  dos  dimensiones  de  la  imagen. 
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Filtrado  digital 


La  imagen  puede  ser  reconstruida  en  el  dominio  es¬ 
pacial  a  partir  de  su  transformada  de  Fourier  apli¬ 
cando  la  transformada  inversa: 


&  i=0  j=0 

(25) 

La  expresión  (22)  también  suele  escribirse 
guíente  modo: 

del  si- 

K- 1  K- 1 

<K(r.s)  =  ZÍ>(¡j)W'“,‘ 

(26) 

i=0  j=0 

donde  W  =  ei2n/K . 

(27) 

La  aplicación  de  la  transformada  de  Fourier  permite 
pasar  la  imagen  del  dominio  espacial  al  de  la  frecuen¬ 
cia.  De  forma  similar  a  lo  que  sucede  con  las  series 
de  Fourier,  también  con  la  transformada  cabe  hablar 
de  espectro  de  magnitud  y  espectro  de  fase  de  la  fun¬ 
ción.  El  sistema  de  coordenadas  en  una  imagen  es  el 
dominio  espacial  de  la  imagen  e  indica  la  posición  de 
cada  celda.  En  el  dominio  de  la  frecuencia,  las  coor¬ 
denadas  ros  representan  la  frecuencia  con  que  un 
determinado  ND  aparece  en  la  imagen. 

El  concepto  de  frecuencia  espacial  de  una  imagen  es¬ 
tática  es  análogo  al  de  frecuencia  temporal  en  el  ám¬ 
bito  del  procesamiento  de  señales.  Una  señal  de  au- 
dio  está  constituida  por  una  composición  más  o  me¬ 
nos  extensa  de  funciones  periódicas.  La  frecuencia 
temporal  indica  la  velocidad  de  cambio  en  los  valores 
instantáneos  de  la  señal.  En  la  imagen,  la  función  a 
estudiar  es  la  sucesión  de  niveles  digitales  correspon¬ 
dientes  a  las  celdas  de  una  misma  línea  o  de  una 
misma  columna.  La  frecuencia  espacial  indica  la  fre¬ 
cuencia  de  alteración  o  cambio  de  los  niveles  digita¬ 
les.  Imágenes  muy  abigarradas  manifiestan  una  alta 


frecuencia  espacial.  En  cambio,  imágenes  muy  unifor¬ 
mes  y  monótonas,  están  caracterizadas  por  su  baja 
frecuencia  espacial.  La  frecuencia  espacial  en  cada  di¬ 
rección  habitualmente  presenta  valores  distintos. 

La  transformada  de  la  imagen  en  dominio  de  la  fre¬ 
cuencia  produce  una  nueva  imagen  <E>(r,s)  en  la  que 
las  coordenadas  de  sus  celdas  representan  la  frecuen¬ 
cia  espacial  en  ambas  direcciones.  La  celda  0(0,0) 
muestra  un  nivel  digital  igual  al  promedio  de  lumi- 
nancia  de  la  imagen,  único  componente  del  espectro 
con  frecuencia  espacial  cero  en  ambas  direcciones. 
Las  demás  celdas  representan  los  componentes  de  la 
imagen  con  frecuencias  que  van  aumentando  en  un 
factor  1/K. 

La  transformada  de  Fourier  tiene  la  propiedad  de  pe¬ 
riodicidad  por  la  cual  el  espectro  de  amplitud  de  la 
imagen  es  simétrico  en  el  entorno  de  K/2  y  el  espec¬ 
tro  de  fase,  antisimétrico.  De  ello  se  deduce  que  la 
imagen  transformada  no  presentará  más  información 
por  encima  de  K/2,  ni  horizontal  ni  verticalmente. 
Por  esto,  se  acostumbra  a  colocar  la  celda  <E>(0,0)  en 
el  centro  de  imagen,  de  tal  modo  que  la  imagen  fre- 
cuencial  resulta  ser  simétrica  respecto  de  su  centro 

Teorema  de  la  convolución 

El  teorema  de  la  convolución  en  el  espacio  establece 
que  la  transformada  de  Fourier  de  la  convolución  en¬ 
tre  dos  funciones  en  el  dominio  del  espacio  es  igual 
al  producto  de  las  transformadas  de  cada  una  de 
ellas.  Enunciándolo  para  una  función  unidimensional 
en  el  dominio  temporal: 

5[/1(t)®/2(t)]  =  F1(co).F2((o)  (28) 

donde  F,  (o>)  =  . ?[/,(*)]  y  F2  (co)  =  -h[/2  (r)]  y  5  es 


FFT5HIFT  3  4 

»  - 

>s.  j  /  :  ;  1 

4  \  ¡  /  3 


Simetría  central  de  la  transformada  de  Fourier  y  aplica¬ 
ción  a  una  imagen  real.  El  origen  de  coordenadas  fre- 
cuenciales  se  sitúa  en  el  centro  de  la  transformada. 


Imagen  real  y  su  transformada  de  Fourier,  con  el  origen 
de  coordenadas  frecuenciales  en  el  centro  del  espectro. 


Carlos  Pinilla 
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la  transformada  de  Fourier. 

De  forma  análoga  se  puede  expresar  el  teorema  de  la 
convolución  en  el  de  la  frecuencia: 

ÍF1  [F,  (a)  0  F2  (a)]  =  271  fx  (t)  f2  (t)  (29) 

siendo  ahora: 


Filtros  de  paso  bajo  y  de  paso  alto  para  aplicar  a 
una  imagen  en  el  dominio  de  la  frecuencia  me¬ 
diante  su  producto. 


F1  [í;  w]  =  /,  (r)  y  5J  [F2  (co)]  =  /,  (») 

La  transformada  inversa  de  la  convolución  de  dos 
funciones  es  igual  al  producto  de  sus  transformadas 
inversas,  afectadas  de  una  constante. 

De  una  forma  sencilla  y  aproximada  se  puede  decir 
que  lo  que  en  el  dominio  del  tiempo  es  una  convolu¬ 
ción,  en  el  de  la  frecuencia  equivale  a  un  producto,  y 
que  un  producto  en  el  dominio  espacial  equivale  a 
una  convolución  en  el  de  la  frecuencia. 


Plantilla  de  filtro  de  paso  bajo  sobre  la  imagen 
frecuencial  y  el  resultado  del  filtrado  devuelto  al 
dominio  espacial  mediante  la  transformada  in¬ 
versa  de  Fourier. 


Del  teorema  de  la  convolución  deriva  una  propiedad 
muy  interesante  de  la  transformada  de  Fourier  apli¬ 
cable  al  filtrado.  El  filtrado  digital,  tal  como  se  ha  ex¬ 
plicado  arriba,  consiste  en  una  convolución  entre  la 
imagen  y  la  función  filtro.  Si  se  pasan  al  dominio  de 
la  frecuencia  tanto  la  imagen  como  el  filtro,  pueden 
ambas  transformadas  multiplicarse  simplemente  para 
realizar  el  filtrado.  Solamente  quedará  para  finalizar 
el  proceso  aplicar  la  transformada  inversa  de  Fourier 
al  producto  para  volver  al  dominio  espacial  y  obtener 
la  imagen  filtrada.  Este  proceso  se  denomina  filtrado 
en  el  dominio  de  la  frecuencia. 


Plantilla  de  filtro  de  corte  de  banda  sobre  la  ima¬ 
gen  frecuencial  y  el  resultado  del  filtrado  de¬ 
vuelto  al  dominio  espacial  mediante  la  transfor¬ 
mada  inversa  de  Fourier. 


La  convolución  z'{i,j)  =  w(m,n)®z(i,j)  propia 

del  filtrado  en  el  dominio  del  espacio  puede  susti¬ 
tuirse  por  el  siguiente  producto  en  el  dominio  de  la 
frecuencia: 


Z'(r,s)  =  W(p,oyz(r,s) 


(30) 


donde  Z'  =  5(z');  W  =  3(w);  Z  =  5(z). 

El  filtro  en  el  dominio  de  la  frecuencia  W  tiene  un 
aspecto  peculiar  y  en  la  mayoría  de  los  casos  viene  a 
ser  un  círculo  central  de  ND=1  en  un  entorno  de 
ND=0,  si  se  trata  de  un  filtro  de  paso  bajo,  o  bien  el 
caso  contrario:  un  círculo  de  ND=0  rodeado  de  ND=1 
hasta  completar  la  imagen.  Como  quiera  que  su  apli¬ 
cación  en  el  dominio  de  la  frecuencia  es  mediante  el 
producto,  cuando  se  multiplique  esta  función  por  la 
imagen  en  el  dominio  de  la  frecuencia  solamente  se 
extraerán  los  valores  centrales  en  el  primer  caso,  o 
solamente  los  periféricos  en  el  segundo.  Dicho  en 
otros  términos,  se  mantienen  las  bajas  frecuencias, 
que  como  se  ha  dicho  se  sitúan  próximas  al  centro 
de  la  imagen  frecuencial,  o  bien  se  mantienen  sola¬ 
mente  las  altas  frecuencias,  en  el  segundo  caso.  De 
ahí  la  denominación  de  filtros  de  paso  bajo  o  de  paso 
alto. 

El  filtrado  en  el  dominio  de  la  frecuencia  es  más  po¬ 
deroso  y  versátil  que  el  filtrado  espacial,  ya  que  per¬ 
mite,  entre  otras  cosas,  dejar  ajenos  al  filtrado  deter¬ 
minados  intervalos  de  frecuencia  espacial  que  con¬ 
vengan  ser  conservados  como  en  la  imagen  original 
Para  evitar  que  en  el  proceso  de  filtrado  de  paso  alto 
se  pierdan  las  características  geométricas  y  radiomé- 
tricas  básicas  de  la  imagen  en  favor  de  los  elementos 
de  sólo  baja  o  alta  frecuencia,  es  posible  idear  filtros 
de  paso  de  banda  o  de  corte  de  banda,  respectiva¬ 
mente.  Estos  últimos  tienen  la  misión  de  extraer  los 
componentes  de  alta  frecuencia  sin  perder  los  valores 
próximos  al  centro  de  la  imagen  frecuencial.  De  este 
modo  se  mantiene  la  imagen  y  se  refuerzan  los  rasgos 
lineales  de  un  modo  similar  a  como  actúa  el  filtro 
laplaciano. 
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Filtrado  digital 


Se  puede  así  mismo  adoptar  para  la  función  de  fil¬ 
trado  frecuencial  valores  intermedios  entre  0  y  1,  de 
modo  que  el  efecto  del  filtro  no  sea  tan  drástico.  Ade¬ 
más,  es  posible  diseñar  filtros  de  muesca  (notch  fil- 
ters )  que  permiten  omitir  determinados  intervalos  de 
frecuencia  en  ciertas  direcciones  espaciales  y  no  en 
otras.  Esto  se  consigue  aplicando  máscaras  de  fil¬ 
trado  frecuencial  de  más  de  una  ventana  en  posición 
no  centrada  con  respecto  al  espectro  de  amplitud  a 
filtrar,  pero  colocadas  simétricamente  con  respecto  al 
centro  del  espectro.  El  resultado  al  deshacer  la  trans¬ 
formación  al  dominio  espacial  es  un  filtrado  direccio- 
nal  pero  sólo  de  los  niveles  digitales  que  muestran  la 
frecuencia  de  cambio  correspondiente  a  la  situación 
de  las  ventanas  muesca  del  filtro. 

Una  vez  realizada  la  operación  de  filtrado  en  el  domi¬ 
nio  de  la  frecuencia,  el  espectro  resultante  debe  ser 
transformado  en  imagen  espacial  mediante  la  aplica¬ 
ción  de  la  transformada  inversa.  En  total,  el  conjunto 
de  las  etapas  de  paso  del  dominio  espacial  al  de  la 
frecuencia,  la  propia  operación  de  filtrado  y  la  rever¬ 
sión  final  al  dominio  espacial  consumen 


Restauración  de  imágenes 

La  adquisición  de  la  imagen  comporta  una  serie  de 
procesos  que  contribuyen  a  la  degradación  de  la  ima¬ 
gen  real.  El  factor  más  importante  suele  ser  la  propia 
atmósfera,  cuyos  aerosoles  en  suspensión  producen 
un  efecto  dispersivo  que  emborrona  la  nitidez  de  la 
imagen.  A  ello  se  suma  la  propia  óptica  del  sensor 
con  sus  sistemas  de  lentes  y  espejos  que  también  dis¬ 
torsionan  la  señal.  Considerando  el  conjunto  de  ele¬ 
mentos  que  degradan  la  señal  como  un  sistema,  la 
transferencia  desde  el  espacio  objeto  al  especio  ima¬ 
gen  vendrá  determinada  por  la  función  del  sistema. 


Ofrec  =2  K2  log2  K  +  K2  multiplicaciones.  En  contra¬ 
posición,  si  el  filtrado  es  realizado  en  el  dominio  es¬ 
pacial  utilizando  una  matriz  deslizante  de  CxL  (co- 
lumnasxlíneas),  el  número  de  operaciones  implicadas 
son  Oesp  =  CLK2.  En  consecuencia,  la  razón  entre  las 
operaciones  necesarias  para  el  filtrado  en  uno  y  otro 
dominio  es:  R  =  Oesp/Ofrec  =  CL¡[2log2  K  +  l).  Para 

una  imagen  de  512x512  celdas  y  una  matriz  desli¬ 
zante  de  filtrado  de  3x3,  el  número  de  operaciones 
en  el  domino  espacial  es  de  2  359  296.  La  misma  ope¬ 
ración  en  el  dominio  de  la  frecuencia  exige  4  980  736. 
La  tendencia  se  afianza  a  medida  que  la  imagen  crece 
de  tamaño.  Sin  embargo,  se  invierte  cuando  se  utili¬ 
zan  plantillas  de  filtrado  de  mayor  número  de  elemen¬ 
tos.  Así  pues,  para  la  misma  escena  anterior,  si  el  fil¬ 
trado  espacial  se  realiza  con  una  matriz  de  4x4  ele¬ 
mentos,  el  número  de  operaciones  resulta  ser  de 
4 194  304,  mientras  que  en  el  dominio  de  la  frecuen¬ 
cia  el  número  de  operaciones  sigue  siendo  el  mismo. 
Si  el  filtro  se  efectúa  con  una  plantilla  de  5x5,  el  fil¬ 
trado  espacial  consume  ya  6  553  600  operaciones. 


De  izquierda  a  derecha  y  de  arriba  hacia  abajo:  imagen  origi¬ 
nal,  su  espectro  en  el  dominio  de  la  frecuencia,  filtro  de 
muesca  y  resultado  del  filtrado  en  el  dominio  espacial.  Obsér¬ 
vese  cómo  se  han  borrado  las  pistas  de  rodadura  y  aterrizaje 
del  aeropuerto,  que  son  los  únicos  elementos  geométricos 
perpendiculares  a  la  orientación  de  las  muescas  del  filtro. 


En  el  capítulo  6  se  introdujeron  los  conceptos  de  fun¬ 
ción  de  dispersión  puntual  ( PSF)  y  del  impulso  uni¬ 
tario.  Pues  bien,  el  patrón  de  degradación  de  la  señal, 
como  ya  se  vio,  viene  establecido  por  la  PSF,  función 
que  describe  la  forma  en  la  que  un  impulso  unitario 
es  transformado  en  su  imagen  final. 

Una  propiedad  del  impulso  unitario  no  vista  hasta 
ahora  es  que  su  transformada  de  Fourier  es  igual  a 
una  constante  de  valor  unidad: 
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(31) 


5[6(í)]=  J  6(í)  e~ia,dt  =  \ 


y  de  un  modo  análogo,  la  transformada  de  Fourier  de 
una  constante  es  un  impulso: 


á(c)  =  |  ce  ,mdt  =  27ic5((o)  (32) 


Ambas  serán  de  utilidad  en  la  restauración  de  imáge¬ 
nes. 

De  no  existir  atmósfera  dispersiva  ni  sistemas  ópticos 
que  introdujeran  distorsión,  la  respuesta  del  sistema 
ante  el  impulso  unitario  debería  ser  otro  impulso  uni¬ 
tario.  Desafortunadamente,  el  sistema  degrada  la  se¬ 
ñal  y  convierte  el  impulso  unitario  en  una  PSF,  que 
es  la  función  dispersiva  h(x,y).  Si  se  conoce  con  pre¬ 
cisión  la  PSF  de  una  imagen,  su  convolución  por  la 
imagen  real  debe  proporcionar  exactamente  la  ima¬ 
gen  finalmente  adquirida.  Siendo  esto  así,  cabe  plan¬ 
tearse  si  fuera  posible  aplicar  una  suerte  de  función 
inversa  que,  aplicada  a  la  imagen  adquirida  pudiese 
proporcionar  la  imagen  real.  Este  es  el  campo  de  la 
restauración  de  imágenes. 

Existen  diversos  métodos  de  restauración  de  imáge¬ 
nes,  de  los  cuales,  unos  operan  en  el  dominio  del  es¬ 
pacio,  otros  netamente  en  el  de  la  frecuencia  y,  final¬ 
mente,  algunos  que  aprovechan  el  dominio  de  la  fre¬ 
cuencia  para  deducir  funciones  restauradoras  que  se 
aplican  sobre  la  imagen  en  el  dominio  del  espacio,  de 
ahí  que  también  se  les  llame  funciones  de  deconvolu- 
ción.  Un  procedimiento  relativamente  sencillo  es  el 
de  deconvolución  ciega  (blind  deconvolution)  que  se 
explica  a  continuación. 

Se  parte  de  la  siguiente  premisa:  si  se  conociera  la 


Restauración  de  imágenes 


10-13 


función  de  deconvolución  g(x,y )  y  la  PSF  del  sistema 
formador  h(x,y)  debería  verificarse: 

g(x,y)®h(x,y)  =  8(x,y)  (33) 

es  decir,  la  convolución  de  la  función  restauradora, 
por  la  PSF  del  sistema,  por  su  propia  naturaleza, 
debe  convertir  esa  PSF  en  un  impulso  unitario, 
puesto  que  su  acción  debe  ser  contrarrestar  el  efecto 
dispersivo  del  sistema.  Esa  misma  expresión,  en  el 
dominio  de  la  frecuencia  será: 

FFT  [g  (x,  y )]  •  FFT  [h (x,  y)]  =  FFT  [ó (x,  y)]  (34) 

siendo  FFT  la  transformada  rápida  de  Fourier,  un  al¬ 
goritmo  para  realizar  la  transformada  discreta  de 
Fourier  sobre  la  imagen  digital.  Escrito  de  otra  ma¬ 
nera: 

G(r,í)-//(r,í)  =  [l]  (35) 

en  donde: 


'1 

1 

1 

...  f 

1 

1 

1 

...  1 

FFT[S(x,q]  =  [l]  = 

1 

1 

1 

...  1 

1 

1 

1 

...  1 

Aplicando  la  propiedad  (29)  a  un  impulso  bidimensio- 
nal.  La  constante  unidad  es  ahora  una  imagen  con- 
tiua  con  todos  sus  ND=1. 

La  función  de  deconvolución  puede  obtenerse  me¬ 
diante  la  reversión  al  dominio  del  espacio  de  la  ex¬ 
presión  (43): 
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Filtrado  digital 


*(*’^FFr'tíy  <3?) 

Para  comprobar  que  la  aplicación  de  g(x,y)  da  los 

resultados  perseguidos,  considérese  que  en  el  pro¬ 
ceso  de  adquisición  de  la  imagen  se  verifica: 

h(x,y)®l(x,y)  =  f(x,y)  (38) 

en  donde  I{x,y)  es  la  imagen  real,  f[x,y)  es  la  imagen 
degradada  y  h(x,y )  es  la  PSF  del  sistema. 

Operando  con  la  función  de  deconvolución  en  los  dos 
miembros: 

g(x,y)®h(x,y)®I  (x,  y)  =  g  (x,  y)®f  (x,  y)  (39) 
o  sea: 

8(x,y)®l(x,y)  =  g(x,y)®f(x,y)  (40) 

Y  como  la  convolución  de  un  impulso  con  una  ima¬ 
gen  da  como  resultado  la  misma  imagen: 

/(^>7)  =  §(^>7)®/(^>7)  (41) 

Expresión  que  viene  a  demostrar  que  la  imagen 
puede  ser  restaurada  mediante  la  convolución  de  la 
función  de  deconvolución  con  la  propia  imagen. 


En  cuanto  a  la  forma  de  obtener  la  PSF  de  la  imagen, 
también  existen  diversos  procedimientos.  Uno  relati¬ 
vamente  sencillo  consiste  en  buscar  dentro  de  la  es¬ 
cena  áreas  de  saltos  bruscos  de  reflectancia,  esto  es, 
zonas  contiguas  de  alto  contraste,  como  pueden  ser 
las  fronteras  agua-tierra  en  el  visible.  Un  perfil  radio- 
métrico  o  transecto  perpendicular  a  la  frontera  de 
ambas  zonas  en  ausencia  de  efecto  dispersivo  debería 
parecerse  a  una  función  escalonada.  Dado  que  la  PSF 
degrada  la  señal  por  su  efecto  dispersivo,  el  escalón 
aparecerá  suavizado.  Las  últimas  celdas  de  la  zona 
oscura  empezarán  a  clarearse,  en  tanto  que  las  pri¬ 
meras  de  la  zona  clara  no  serán  tan  brillantes  como 
debieran.  La  derivada  de  esa  curva  de  transición  es 
precisamente  la  PSF  buscada. 


Celda 


Arriba:  transectos  sobre  una  presa,  zona  de 
alto  contraste  en  su  cara  aguas  arriba. 
Abajo,  la  PSF  obtenida  al  derivar  el  prome¬ 
dio  de  los  perfiles  radiométricos  obtenidos. 
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